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Time reversal protected three dimensional (3D) topological paramagnets are magnetic analogs of 
the celebrated 3D topological insulators. Such paramagnets have a bulk gap, no exotic bulk exci¬ 
tations, but non-trivial surface states protected by symmetry. We propose that frustrated spin-1 
quantum magnets are a natural setting for realising such states in 3D. We describe a physical picture 
of the ground state wavefunction for such a spin-1 topological paramagnet in terms of loops of fluc¬ 
tuating Haldane chains with non-trivial linking phases. We illustrate some aspects of such loop gases 
with simple exactly solvable models. We also show how 3D topological paramagnets can be very 
naturally accessed within a slave particle description of a spin-1 magnet. Specifically we construct 
slave particle mean field states which are naturally driven into the topological paramagnet upon 
including fluctuations. We propose bulk projected wave functions for the topological paramagnet 
based on this slave particle description. An alternate slave particle construction leads to a stable 
U(l) quantum spin liquid from which a topological paramagnet may be accessed by condensing the 
emergent magnetic monopole excitation of the spin liquid. 


Frustrated quantum magnets display a rich variety of 
many-body phenomena. Some such magnets show long- 
range magnetic order at low temperature, often selected 
out of a manifold of degenerate classical ground states 
by quantum fluctuations. A very interesting alternative 
possibility — known as quantum paramagnetism — is 
the avoidance of such ordering even at zero tempera¬ 
ture. Quantum paramagnets may be of various types. 
A fascinating and intensely-studied class is the quan¬ 
tum spin liquids: these display many novel phenomena, 
for instance fractionalization of quantum numbers and 
topological order, or gapless excitations that are robust 
despite the absence of broken symmetries [IH3]. 

Recently there has been much progress in understand¬ 
ing a different type of remarkable quantum paramagnet. 
These are phases which have a bulk gap and no fractional 
quantum numbers or topological order. Despite this, 
they have nontrivial surface states that are protected 
by global symmetries. These properties are reminiscent 
of the celebrated electronic topological band insulators. 
Hence they have been called topological paramagnets [ 1 ] . 
Topological paramagnets and topological band insulators 
are both examples of what are known as Symmetry Pro¬ 
tected Topological (SPT) phases [SHZ]. A classic example 
of a topological paramagnet is the Haldane/AKLT spin-1 
chain: though this has a bulk gap and no bulk fraction¬ 
alization, it has dangling spin-1/2 moments at the edge 
which are protected by symmetry, for instance time re¬ 
versal. In the last few years tremendous progress has 
been made in understanding such SPT phases and their 
physical properties in diverse dimensions (for reviews, see 
Refs. IHlandini). 

The main focus of the present paper is on three- 
dimensional topological paramagnets that are protected 
by time reversal (we also briefly discuss topological para¬ 
magnets protected by other symmetries, notably conser¬ 
vation of at least one spin component). These are inter¬ 
esting for a number of reasons. First, time reversal is a 
robust symmetry of typical physical spin Hamiltonians. 


In ID the familiar Haldane/AKLT chain is the only time 
reversal protected topological paramagnet while in 2D 
there are no time reversal protected topological param¬ 
agnets. In 3D however there are three distinct non-trivial 
phases [U [101 HI] (corresponding to a classification by 
the group Z|). Second, regarded as an electronic insu¬ 
lator, unlike the ID Haldane chain [12], these 3D topo¬ 
logical paramagnets survive as distinct interacting SPT 
insulators m- The properties and experimental finger¬ 
prints of such topological paramagnets were described in 
Refs. jHHni dH US] • However there is currently very lit¬ 
tle understanding of where such phases might actually be 
found. In this paper we propose that frustrated spin-1 
Mott insulators may be good places to look for an exam¬ 
ple of such phases. 

Already in the familiar ID example it is the spin-1 anti¬ 
ferromagnetic chain, rather than the spin-1/2 chain, that 
naturally becomes a topological paramagnet. In 3D for 
one of the topological paramagnets we provide a phys¬ 
ical picture and a parton construction which are both 
very natural for the spin-1 case. We hope that our ob¬ 
servations inspire experimental and numerical studies of 
frustrated spin-1 quantum magnetism in the future. To¬ 
wards the end of the paper we remark on materials that 
may form such interesting frustrated magnets. 

The three 3D topological paramagnets that are pro¬ 
tected by time reversal symmetry alone n dni E] all 
allow for a gapped surface with Z 2 topological order {i.e. 
a gapped surface Z 2 quantum spin liquid) even though 
the bulk itself is not topologically ordered. The prop¬ 
erties of this surface theory give a useful way to label 
the bulk phases. The surface has gapped quasiparticle 
excitations — labelled ‘e’ and ‘to’ — which are mutual 
semions. These may be thought of as the electric charge 
and magnetic flux of a deconfined Z 2 gauge theory (like 
the vertex and plaquette defects of Kitaev’s toric code 
[14]). At the SPT surfaces these particles have properties 
— self-statistics or time reversal transformation proper¬ 
ties — that are impossible in a strictly 2D system, and 
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which encode the topology of the bulk wavefunction. The 
three nontrivial bulk states are denoted: 

eTmT, efTmfT, efmf. 

In the first and second, the surface e and m excitations 
are each Kramers doublets under time reversal, denoted 
by T. In the second and third they are fermions (/), 
while in the first they are bosons. This paper focuses 
primarily on the ‘eTmT’ state. 

We begin by explaining a physical picture of a suit¬ 
able ground state wave function for the eTmT topologi¬ 
cal paramagnet. This is most easily visualized on a dia¬ 
mond lattice. We first close-pack each interpenetrating 
fee sublattice of the diamond lattice with closed loops. 
On each loop we place all the spin-1 moments (located at 
the diamond sites) in the ground state of the ID AKLT 
chain. We then superpose all such loop configurations 
with a crucial (—1) sign factor whenever loops from the 
two different fee sublattices link. We argue that this con¬ 
struction yields the topological paramagnet. 

To understand the topological properties of such a 
wave function we describe a simple exactly solvable loop 
gas Hamiltonian m — equivalent to two coupled Ising 
gauge theories — that clarifies the role of the ‘(_i)i'nking) 
sign structure. In this solvable model the loops do not 
have AKLT cores but there are two species of loops on 
different sublattices with the mutual (—1) linking sign. 
It demonstrates very simply how this sign leads to a state 
without intrinsic topological order. (This loop gas is not 
in the eTmT state, because of the absence of AKLT 
cores, but we show it to be nontrivial in a different sense.) 

Next we use the two-orbital fermionic parton repre¬ 
sentation developed for spin-1 magnets m to construct 
possible ground states. When the fermionic partons have 
the mean-field dispersion of a certain topological super¬ 
conductor, we show that the gauge fluctuations associ¬ 
ated with the parton description convert the system into 
a topological paramagnet. In this construction the mean 
field state is unstable toward confinement by gauge fluc¬ 
tuations, as a result of a continuous nonabelian gauge 
symmetry. Despite this the bulk gap survives, leaving 
behind a non-trivial surface that we are able to iden¬ 
tify as that of the eTmT topological paramagnet. As a 
warm up exercise to illustrate some of the ideas of this 
3D construction, we also describe how to access the ID 
Haldane phase by confining a topological superconduc¬ 
tor of parton fermions. The 3D construction naturally 
suggests alternative bulk wave functions for topological 
paramagnets, in the form of Gutzwiller-projected topo¬ 
logical superconductors. This may be fruitful for future 
numerical work on the energetics of microscopic models. 

This parton construction also gives access to other SPT 
states for quantum magnets in 3D. For instance we show 
how to naturally obtain an SPT paramagnet (dubbed 
eCmT in Ref. [TOl) protected by U(l) x Zl’, where the 
U(l) describes rotation about one spin axis, say Sz, and 
2-2 is time reversal. 

Finally we show how to access a bulk U(I) quantum 


spin liquid with non-trivial implementation of time rever¬ 
sal symmetry. Interestingly simply condensing the mag¬ 
netic monopole of this U(I) spin liquid leads to an SPT 
state dubbed eCTmT in the presence of both spin rota¬ 
tion and time reversal symmetries. If only time reversal 
is present this becomes the eTmT state. 


I. LOOP GAS STATES 

In this section we describe a loop gas wavefunction 
that is naturally adapted to spin one magnets and gives 
an intuitive picture for the eTmT state. The wavefunc¬ 
tion is a superposition of loop configurations, with each 
loop representing an AKLT state m for the spins ly¬ 
ing on it. A given configuration enters the superposition 
with a sign factor determined by its topology: specifi¬ 
cally, the loops come in two species, A and B (one as¬ 
sociated with each sublattice of the bipartite diamond 
lattice) and the sign depends on the linking number of A 
loops with B loops. This geometrical picture makes the 
relationship between the bulk wavefunction and the sur¬ 
face excitations particularly simple. The surface e and 
m excitations are endpoints of the two species of AKLT 
chains, and are Kramers doublets since an AKLT chain 
has dangling spin-1/2s at its ends. 

In Sec. [H] we describe a similar wavefunction for ‘pure 
loops’, i.e. loops that do not carry an internal AKLT 
structure. This may be regarded as a state of two coupled 
Ising gauge theories. It is not in the eTmT phase, but it 
illustrates the basic features of the loop gases in a simple 
model with an exactly solvable Hamiltonian. This ‘pure 
loop’ model is also interesting in its own right: when open 
strands (as opposed to closed loops) are banished from 
the Hilbert space, i.e. when charge is absent, it is in a 
nontrivial phase despite the absence of topological order. 
Therefore it may be viewed as a ‘constraint-protected’ 
state. It would be interesting to relate this to the recent 
ideas of Ref. [18]. We note that the models discussed 
in Ref. [Tn| are also believed to be separated from the 
trivial phase by a phase transition, despite the absence 
of topological order. 

The wavefunctions discussed here are in a similar spirit 
to the Walker Wang models, which are formulated in 
terms of string nets with a nontrivial sign structure, 
and show bulk confinement and surface topological or¬ 
der [TTl 1201 mi- Constructions of SPTs using Walker 
Wang models were given in Refs. mm- 2D ‘symmetry- 
enriched’ topological states and SPT states [ 20 ] 

have also been constructed by attaching AKLT chains to 
loop-like degrees of freedom (see also [ 17 ]). 

A. Fluctuating AKLT chains 

The diamond lattice is made up of two fee sublattices, 
A and B. If Ca is a configuration of fully packed loops 
on A (with every A site visited by exactly one loop), we 
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define \Ca) to be a product of AKLT states \C) for each 
of the loops C in Ca, 

\ca) = n i^) • (1) 

C^Ca 




Similarly \Cb) is the state corresponding to a loop config¬ 
uration Cb on B. To define the AKLT states |£) fully we 
must choose an orientation for the fee links, as discussed 
below (Sec. IB|. 

Let X (Ca , Cb ) be the mutual linking number of the two 
species of loops. Since the loops are unoriented, this is de¬ 
fined modulo two: X{Ca,Cb) = 0,1. A schematic wave- 
function for the eTmT phase may be written in terms of 
X{Ca,Cb). 


1‘J’) = E |Ca) \Cb) ■ (2) 

Ca,Ca 

For concreteness, we take periodic boundary conditions. 
The sums over Ca and Cb are then each restricted to loop 
configurations with an even number of strands winding 
around the 3D torus in each direction, for reasons dis¬ 
cussed below. This global constraint, together with the 
geometrical fact that the links of A never intersect those 
of B, ensures that X{Ca,Cb) is well defined. 

The entanglement between the two sublattices in Eq. 
is entirely due to the sign factor. First consider what hap¬ 
pens in the absence of this sign factor. Each sublattice 
then hosts a superposition of loop configurations with 
positive amplitude, e.g. J2ca I^^)- analogy with the 
usual picture of deconfined Z 2 gauge theory as a super¬ 
position of electric flux loop configurations [5S] , we would 
expect such a state to show Z 2 topological order. (It is 
a 3D version of the ‘resonating AKLT’ states studied in 
2D [251 - 55] . 1 The endpoint of an open AKLT chain is 
the deconfined Z 2 charge in this state. Associated with 
the topological order is ground state degeneracy — dif¬ 
ferent ground states are distinguished by the parity of 
the winding number in each spatial direction. 



FIG. 1. Two species of AKLT loops, one on each sublattice 
of the diamond lattice. The loops live on the links of the fee 
sublattices, i.e. on next-nearest-neighbour bonds of diamond. 


FIG. 2. For appropriate boundary conditions, endpoints of A 
and B chains (red and blue respectively) give surface excita¬ 
tions with mutual semionic statistics. Braiding the anyons on 
the surface (first arrow) changes the sign of the wavefunction, 
for consistency with the rule that configurations related by 
passing an A strand through a B strand in the bulk (second 
arrow) appear in the wavefunction with opposite sign. 


In contrast, |$) is not expected to show topological 
order, despite the proliferation of long loops in Eq. 
Instead it describes a phase in which the endpoints of 
open chains are confined in the bulk. Furthermore there 
is no ground state degeneracy: states with odd winding 
numbers are not ground states (i.e. are not locally indis¬ 
tinguishable from 1$)). 

More detailed discussion of this is deferred for the solv¬ 
able model of Sec. [n] but the basic idea is the following. 
While the amplitude (—1)^(^-4 Tb) depends on the global 
topology of the loop configurations, it amounts to the 
simple loeal rule that the amplitude changes sign if an A 
strand is passed through a B strand. It is useful to imag¬ 
ine a hypothetical parent Hamiltonian that imposes this 
sign rule. But the sign rule cannot be consistently im¬ 
posed if the wavefunction includes open strands or config¬ 
urations with odd winding numbers (see below). Similar 
phenomena occur in the confined Walker-Wang models 

[iiiiiniiii]. 

However, open endpoints are deconfined at the bound¬ 
ary, for appropriate boundary conditions. The minus 
sign associated with passing an A strand through a B 
strand in the bulk means that the endpoints are mutual 
semions [29] — see Fig. They are also Kramers dou¬ 
blets. These surface properties are the defining features 
of the eTmT state. The wavefunction |$) has more sym¬ 
metry than simply time reversal (e.g. separate spin rota¬ 
tion symmetries for each sublattice) but if it is indeed in 
the eTmT phase then these symmetries could be weakly 
broken without leaving the phase. 


B. Further details on fluctuating AKLT state 

To write the AKLT-based state explicitly it is conve¬ 
nient to represent the spin-one at each site i in terms of 
auxiliary spin-1/2 bosons [T7|[27]. If the boson creation 
operators are bj^ (a then = \b\^aapbip. The 

occupation number h\^bia is equal to two to ensure spin 
one at each site. The AKLT state |£) is then created 
by acting on the boson vacuum with operators that 
create singlet pairs on the links of the loop, which we 
normalize as 3}^ = ~ operator 













4 



FIG. 3. Left: Loops on interpenetrating cubic lattices A and 
B. The state |’I') is a snperposition of such configurations 
with signs determined by linking of A and B loops. Right: 
the prodnct of Panli matrices defining the flip term 7^ on a 
plaquette (see Eq. [^. 


is antisymmetric in (i,j), so to define |$) we must fix 
an orientation for the links of each fee sublattice. (The 
fee lattice has four sublattices, a, b, c, d, so for example 
we could orient the links from a ^ a ^ c, a ^ d, 
b ^ c ^ d ^ b, with the orientations on each sublattice 
related by inversion symmetry.) Then for each sublattice 

\c)= n 4 i^ac), (3) 

(b>ec 

where i is the site at the tail of the oriented link {ij). 
These states satisfy (C|C) = nioops(l + (“1)^/3^”^): 
where ^ is the length of a given loop un¬ 
it should be noted that that expectation values in the 
state 1$) are nontrivial, in particular because overlaps 
{C\C') for distinct C, C are nonzero. So while it is plausi¬ 
ble that 1$) is in the eTmT phase, this cannot be estab¬ 
lished purely analytically. For example, the state could in 
principle break spatial or spin rotation symmetry sponta¬ 
neously. A cautionary example is given by the uniform- 
amplitude resonating valence bond state for spin-l/2s on 
the cubic lattice: this has weak Neel order m, despite 
being a superposition of singlet configurations which in¬ 
dividually have trivial spin correlations. In the present 
model, the entanglement between sublattices supresses 
off-diagonal elements of the reduced density matrix when 
written in the AKLT-chain basis |3T]. Together with the 
non-bipartiteness of the fee lattice, this makes spin order 
seem less likely. But since |<I>) is intended to illustrate the 
topological structure of the phase, and not as a ground 
state of a realistic Hamiltonian, it may not be crucial 
whether it is in the desired phase as written or whether 
further tuning of the amplitudes is required. 


II. ‘PURE LOOP’ STATE 

It is enlightening to look at the simplest model[Tn] that 
captures the (—sign structure. To this end we 
take a system of spin-I/2s on the links of two interpen¬ 
etrating cubic lattices A and B, as shown in Fig. We 
think of a down spin (in the ‘z’ basis) as an occupied 


link, and an up spin as an unoccupied one. The num¬ 
ber of occupied links at each vertex is always even in the 
state we consider, so the configurations of occupied links, 
Ca and Cb, can be decomposed into closed loopsQ We 
refer to Ca and Cs as loop configurations. Other solv¬ 
able loop gas/string net models have been considered in 
Refs, [m HI] , using the Walker Wang construction HD] ■ 
The ‘pure loop’ state analogous to |$) above is (again 
we sum only over loop configurations with even winding 
numbers on each sublattice): 

I'f') = E \Ca) \Cb) . (4) 

Ca,Cb 


We may view Ca and Cs as the electric flux line config¬ 
urations for a pair of coupled Z 2 gauge fields, with one 
Z 2 gauge field living on each cubic lattice. Imposing the 
above sign structure for the two sets of electric flux lines 
is equivalent to binding the electric flux line of each gauge 
field to the magnetic flux line of the other, as will be clear 
shortly. 

It is straightforward to write down a gapped parent 
Hamiltonian Flunking for j'k), using the fact that flipping 
the occupancy of all the links on the plaquette changes 
the linking number X{Ca,Cb) if and only if the link pierc¬ 
ing the plaquette is occupied. Flunking is a sum of terms 
for the plaquettes p of each cubic lattice: 


Flunking — I ^ ^ ^ d' Ap T J ^ ^ d~ Bp | • (b) 

p&A peB 


The operators J-a and Xb flip the occupancy of the links 
on a plaquette, with a sign that depends on whether the 
link piercing it is occupied. Allowing p to denote both a 
plaquette and the link piercing it, and denoting the Pauli 
operators on A and B hy a and r respectively, 

J^Ap = r; n ^Bp = a; n rf. (6) 

lep lep 

These operators all commute, so the Hamiltonian is triv¬ 
ially solvable, jlk) is the unique ground state and min¬ 
imises each term of Flunking since F^|4') = |4>) for each 
plaquette operator. 

The state |'k) contains only closed loops, i.e. it satisfies 

Y[crf = 1 for u G A, = ^ for u G R (7) 

iGv iGv 


where v denotes a vertex and I G v the links touching v. 
Any state satisfying |'I') = jT) for all the plaquette op¬ 
erators must also satisfy these vertex conditions, because 
riiGu riiG-u '''i written as products of Ts. 


^ With a harmless ambiguity when the number of occupied links 
at a vertex exceeds two. 
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FIG. 4. After a basis change, l^t) is a superposition of mem¬ 
brane configurations = — 1 on shaded plaquettes) with red 
loops (where = — 1) glued to membrane boundaries. (The 
red loops are cr-electric lines and the membrane boundaries 
are r-magnetic lines.) 


We may regard Eqs. as the gauge constraints for a 
pair of pure Z 2 gauge theories (the Z 2 versions of W.E = 
0). The two electric fields are given by and and 
live on the links of A and B respectively. The magnetic 
field of each gauge field lives on the links of the opposite 
lattice to its electric field. For example the magnetic field 
of a is given by Oigp ^ where p is a plaquette of A, or 
equivalently a link of B. 

In this language, Hunking simply glues the electric flux 
line of each species to the magnetic flux line of the other. 
The CT~magnetic flux and the r-electric flux are equal 
since Ba = 1, and the cr-electric and r-magnetic fluxes 
are equal via Eb = 

The state |T') is not topologically ordered. Neither is 
it a time-reversal protected SPT: it can be adiabatically 
transformed to a product state without breaking time 
reversal symmetry. However it is protected if impose 
Eqs. as constraints: i.e. if we forbid open strands, as 
opposed to closed loops. In the gauge theory language, 
this means forbidding charge. With this constraint it is 
impossible to reach a trivial state without going through 
a phase transition, as follows from the self-duality of the 
state described in Sec. EAl 

We will explain these features from several points of 
view below. One convenient approach which leads to a 
geometric picture is to switch from the (cr^, r^) basis used 
in Eq. I^to the ((t*,t^) basis. The cr^ configuration is a 
loop configuration on the A lattice, as above. We rep¬ 
resent the configuration by a configuration of mem¬ 
branes made up of plaquettes on the A lattice: = — 1 

represents an occupied plaquette, and = 1 an unoccu¬ 
pied one. 

The ipB terms in "Hunking act on a link of the A lattice 
together with the four plaquettes touching it. = 1 
imposes the rule that the loops are glued to the bound¬ 
aries of the membranes, i.e. to the links where an odd 
number of occupied plaquettes meet. This is the gluing 
of cr-electric flux lines (where tr^ = —1) to r-magnetic 
flux lines (where = mentioned above. 

Let A4 denote a membrane configuration, and \Ai) 


the corresponding state with = — 1 on the occupied 
plaquettes. Let dM. be the loop configuration given by 
the boundaries of the membranes in M.. Then Irk) can 
be written (neglecting an overall constant) 

l^) = E E \Ca)\M). (8) 

Fig.i shows the geometrical interpretation of this state. 
It is a soup of r“ membranes, with loops glued to their 
boundaries. 

Conhnement of string endpoints is easy to see in this 
basis. A pair of vertex excitations at which JI/gj; '^i ~ 
are connected by an open string. Since the boundary of 
M. contains only closed loops, the open string makes it 
impossible to satisfy the gluing of strings to membrane 
boundaries demanded by the Eb terms in Hunking- If the 
separation of the vertex defects is H, there must be at 
least D unsatisfied links, giving a linear confining poten¬ 
tial for such defects. For similar reasons, a configuration 
with an odd number of winding strands in some direc¬ 
tion costs an energy proportional to the spatial extent of 
the system in this direction. By symmetry, this applies 
equally to the strings that are present in the original 
basis. 

We can also understand the confinement of string end¬ 
points algebraically (Refs. [TTl HT] give analogous argu¬ 
ments for bulk confinement and surface topological order 
in the Walker Wang models). The Hamiltonian in Eq. 
is clearly exactly soluble not just for the ground state 
but for all excited states. An ‘elementary’ excitation is 
given by a ‘defect’ in some square plaquette, say on the 
B lattice, with 


Bbp =-I (9) 

while J- = -1-1 on all other plaquettes of either sublattice. 
Such a defect plaquette costs energy 2J. It leads to a 
violation of the closed loop vertex constraint for cr* on 
the two vertices of the A sublattice connected by the 
A-link that penetrates the defect plaquette. Thus the 
excitation we have created has two string end-points on 
nearest neighbor A-sites. To move these string endpoints 
apart by a distance D we must create o{D) such defect 
plaquettes. Consequently the energy cost is also o{D) 
and we have linear confinement of string endpoints. 

In the gauge theory language, the reason for the ab¬ 
sence of deconfined excitations is that the tensionless 
lines in this state are not lines of pure electric flux, but 
rather of electric flux together with magnetic flux of the 
other species. If such lines could end, their endpoints 
would be deconfined excitations. But the Hilbert space 
does not allow for such excitations: a magnetic flux line 
cannot terminate in the bulk (by virtue of its definition 
in terms of e.g. 

Despite the lack of deconfined endpoints in the bulk, 
A and B strings that terminate on a boundary can give 
deconfined e and m particles in a surface Z 2 topologically 
ordered state. To see this, we terminate the system as in 
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Fig.§ including in the Hamiltonian the natural plaque- 
tte and vertex terms at the surface. The surface string 
operators that create pairs of e or pairs of m excitations 
can then be written explicitly (see Fig. [^. They sat¬ 
isfy the same algebra as the string operators in the toric 
code |14j . confirming that e and m are mutual semions 
as expected from the heuristic argument of Fig. 

We can adiabatically transform to a product state 
so long as we allow the intermediate states to violate the 
closed-loop constraints on at least one sublattice. The 
membrane picture gives an obvious way to do this, by 
giving the membranes in A4 a surface tension. If ‘Area’ 
denotes the number of occupied plaquettes in Ad, the 
interpolating state is 

I '^)7 = E E ( 10 ) 


When 7 = 0 this is the initial state, and when 7 —>■ oo 
only the term with zero area survives. This is the state 
with no loops and no membranes, i.e. the product state 
\a^ = 1) \t^ = 1). To get a gapped parent Hamiltonian 
for we modify the plaquette flip term Ta in 'Hiinking 


to Tap = (cosh 7 ) 


-1 


TpU 




. This pre- 


■ (sinh 7 ) Tp 

serves the simple algebraic properties of the plaquette 
terms. 


A. Self-duality of jd') and protection by constraints 

When the interpolating state above is rewritten in the 
original (cr^,T^) basis, it includes configurations with 
open strands, as well as closed loops, on the B lattice. 
What if we impose the constraint that both lattices have 
only closed loops? In this case it is impossible to go from 
I'l') to a trivial state without a phase transition. (We will 
take the reference trivial state to be that with no loops, 
Itrivial) = |cr^ = I) \t^ = I).) 

This follows from a simple duality transformation 
which exchanges the electric flux of each species with the 
magnetic flux of the other species. The duality maps |d>) 



FIG. 5. String operators creating surface excitations. Left: 
acting with a chain of operators on the links of the upper 
layer {A lattice surface) gives a pair of e excitations (i.e. end¬ 
points of bulk A strings). Right: a pair of m excitations (i.e. 
endpoints of B strings) are created by a chain of operators 
(thick green strand) on the lower layer {B surface), together 
with (7^ operators on the corresponding links in the upper 
layer (thick purple links). 



FIG. 6. Under the mapping (111, a cr^ (or r^) operator on a 
link is exchanged with a product of (resp. cr’’’) operators 
on the surrounding links of the other lattice. (Links of one 
lattice can equally be thought of as plaquettes of the other.) 


to itself, but exchanges the trivial state with a topologi¬ 
cally ordered one. Thus there is no adiabatic path from 
Ilk) to the trivial state. If there were, duality would yield 
an adiabatic path from |'k) to the topologically ordered 
state, and this is impossible since Idl) is not topologically 
ordered. 

The duality transformation makes sense for states 
obeying the closed loop constraint. (To be precise, we 
must also impose the global constraint that the loop 
configurations have even winding in each direction.) As 
shown in Fig. its action is: 

( 11 ) 

pel lep 

Here p G I denotes the four plaquettes p surrounding link 
1. We have labelled the as by I for link and the rs by p for 
plaquette, but the duality acts on the two sets of degrees 
of freedom symmetrically. It preserves the locality of any 
Hamiltonian acting in the constrained Hilbert space. 

For completeness, we write the action of the duality 
on states explicitly. Return to the picture of loops -I- 
membranes on the A lattice, i.e. the {a^, t^) basis. One 
may check that any state satisfying the constraints can 
be written as a sum over two loop configurations on the 
same lattice, 

I/) = ^ /(Ca,C:i)|C^)J^),, (12) 

CAfi'^ 


where \C'j^ is defined as the uniform superposition of all 
membrane configurations |Ad)^ with boundary dM. = 
We have added subscripts to the kets as a reminder 
of the degrees of freedom involved. {Ca is the a-electric 
flux configuration, and C\ the r-magnetic flux config¬ 
uration; the fact that the wavefunction depends on M. 
only through dM. is simply a statement of gauge invari¬ 
ance.) The duality then simply exchanges the two kinds 
of loops. 


f{CA,C'^)^f{CM- ( 13 ) 

The flip operators Ba and Fb (Eq. [ 6 ) are clearly in¬ 
variant under the duality in Eq. El and therefore so is 
??ii nking - (We Can also see that jlk) is invariant from 
Eq. 13 and Eq. |^) On the other hand, the trivial Hamil¬ 
tonian 


^trivial — ( J ^ T </ ^ T, 


leA 


leB 


(14) 
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is exchanged with 

T^deconfined — j ^ (15) 

y PGj4 igp pG-B ZGp 

which describes a pair of deconfined Z 2 gauge theories. 
This establishes the claim at the beginning of this sub¬ 
section: while the linking state is invariant, the trivial 
state is exchanged with a topologically ordered state. It 
follows that the linking state is in a distinct phase from 
the trivial state if we do not allow open endpoints in the 


Hilbert space. (We know from Eq. 10 that they are in 


the same phase if we do allow endpoints.) 


B. Heuristic relation between symmetry protection 
of eTmT and closed—loop constraint 

The proposed wave function for the eTmT phase 
has the two loop species ‘stuffed’ with Haldane/AKLT 
chains. The linking sign factor ensures that the ground 
state is not topologically ordered as required for a topo¬ 
logical paramagnet. In particular the open end-points of 
the loops — which now harbor a Kramers doublet — are 
confined. However as described in Sec. ira the surface 
implements time-reversal ‘anamolously’ exactly charac¬ 
teristic of the eTmT state. 

We now briefly consider whether the results in the pre¬ 
vious subsection for the ‘pure loop’ state yield a heuris¬ 
tic ‘bulk’ understanding of why the eTmT state is pro¬ 
tected by time reversal. So let us imagine perturbing the 
schematic eTmT wavefunction of Sec. m and ask why 
we cannot reach a trivial state without a phase transition. 

We make use of the heuristic analogy between the 
AKLT loops of the spin-1 system and the ‘pure loops’ of 
the coupled gauge theory]^ The result for the pure loop 
state then indicates that if we only have closed AKLT 
loops on each sublattice, we cannot get to a trivial state 
without a phase transition. So, we must consider pro¬ 
liferating open strands on at least one sublattice. But 
in the spin-1 system, unlike the pure-loop system, open 
strands introduce bulk spin-1/2 Kramers doublet degrees 
of freedom. (Binding these emergent spin-l/2s into sin¬ 
glets with others on the same sublattice merely heals the 
AKLT chains, taking us back to the original situation 
with separate closed loops on each sublattice.) When 
time reversal is broken, these spin-l/2s are innocuous 
— for example we can gap them out using a magnetic 
field. But it is natural to expect that when time rever¬ 
sal is preserved they prevent us reaching a trivial state 
without closing the gap. 


^ Recall that the singlet basis allows us to represent any spin-zero 
state of the spin-1 system in terms of loops of spin-1/2 singlet 
bonds; these may form minimal length ‘loops’ which backtrack 
on a single link, i.e. spin-1 singlet bonds, or longer AKLT loops. 


However, the above argument is incomplete as it does 
not rule out the possibility of getting to a trivial state 
by proliferating nearby pairs of open strands on opposite 
sublattices. Such a pair gives two spin-l/2s which can be 
bound into a singlet to avoid a gapless degree of freedom. 
In the gauge theory, such pairs correspond to bound pairs 
of electric charges, one from each Z 2 gauge field. The sta¬ 
bility of the eTmT state suggests that the pure loop state 
remains protected even when such double charges are al¬ 
lowed. We note that at the surface these double charges 
correspond to the bound state of the e and m particle (in 
the surface topological order). This is a Kramers singlet 
spin-0 fermion (conventionally denoted e). The surface 
Fermi statistics suggests a potential obstruction to ‘trivi¬ 
alizing’ the bulk by proliferating the double charges. We 
leave an explicit demonstration of this for the future. 


III. BARTON CONSTRUCTIONS 

Though the description of the eTmT topological para¬ 
magnet in terms of a loop gas wave function is physically 
appealing it is desirable to have alternate descriptions 
which enhance our understanding and which may help 
with evaluating the energetic stability of this phase in 
microscopic models. To that end, in this section we pro¬ 
pose explicit parton constructions for some topological 
paramagnets in spin-1 systems. 

Historically the parton approach has provided varia¬ 
tional wave functions and effective field theories both 
for spin liquids [3] and non-fractionalized symmetry¬ 
breaking states |32| . The parton construction inevitably 
introduces a gauge symmetry. It describes a fraction¬ 
alized spin liquid phase whenever it yields an emergent 
deconfined gauge held. To obtain a non-fractionalized 
phase such as conventional antiferromagnet or valence 
bond solid paramagnet, the gauge held should either be 
Higgsed or conhned. 

Recently the parton construction has been used to con¬ 
struct SPT states in two |33l |34l and three [35l |36] di¬ 
mensions. The general idea is to construct a gauge the¬ 
ory (with matter helds) that is conhned, but with cer¬ 
tain non-trivial features surviving in the conhned state 
that make it an SPT state. However, the currently 
known constructions in three dimensions use either Z 2 
or U(l) gauge theories, which do not conhne automati¬ 
cally: strong gauge coupling is needed to reach the con¬ 
hned phase. Furthermore, the constructions using U(l) 
gauge theories [35j require highly nontrivial dynamics of 
the gauge helds to condense composite dyon-like objects. 

In three dimensions, a continuous non-abelian gauge 
symmetry is needed to guarantee conhnement. We pro¬ 
pose two parton constructions in three dimensions with 
SU(2) gauge symmetry, which conhne even if the bare 
gauge coupling is small, giving rise to topological param¬ 
agnets. A similar construction was used previously [34] 
in 2D to describe an SPT phase of a spin-1 magnet pro¬ 
tected by spin SU(2) symmetry and time reversal. We 




also propose a construction with U(l) gauge symmetry, 
which confines at sufficiently strong coupling. Crucially, 
this U(l) construction differs from previous ones in that 
we only condense simple monopoles to confine the gauge 
theory, which can be achieved at strong coupling without 
exotic form of gauge field dynamics. 

The spin-1 operators are re-written using the two- 
orbital fermionic parton representation proposed in 
Ref. [ini 

'^= d (16) 

^ a=l,2 

where a = 1, 2 is the orbital index. As will be discussed 
below, the two-orbital structure is natural for topological 
bands corresponding to topological paramagnets. This 
gives another reason for favoring spin-1 systems. 

The physical spin states are represented in the parton 
description as 

I'S'z = 6) = ^ /it■^2^) 

|S'z = +1) = \S^ = -1) = /i\/selvae). 

where |vac) is the state with no fermions. States in the 
physical spin Hilbert space thus have two fermions at 
each site, J2aa flafaa = 2, and the two fermions form a 
singlet in orbital space: denoting the Pauli matrices in 
orbital space by this is y\ f^^fabfba = 0. 

The representation in Eq. |16| actually has an Sp(4) 
gauge redundancy [T^ which becomes apparent when 
we represent the fermions using Majoranas, / = 

i (? 7 i — ir] 2 ). Here 771^2 are Hermitian operators satisfy¬ 
ing {risi,r]s> j} = 2dss'Sij, where s, s' = 1,2 are the new 
indices associated with the Majoranas and /, J represent 
all other indices (site, spin, orbital). The Majorana rep¬ 
resentation of the spin is 

5=ir7^Er7, S = (pV^a^pV), (17) 

where are Pauli matrices acting on the Majorana 

index. The generators of the gauge symmetry are ten 
anti-symmetric imaginary matrices that commute with 
the physical spin operators: 


r = (18) 


where are Pauli matrices acting on the orbital index. 
The spin in Eq. 17 is invariant under the Sp(4) gauge 
transformation 77 —e“‘^"? 7 . 

The effective field theory associated with the parton 
construction is a gauge theory. The gauge symmetry is 
determined by the mean field band structure of the par- 
tons, and is in general a subgroup of the full Sp(4) group 
due to some generators being Higgsed. The gauge struc¬ 
ture allows symmetry to act projectively on the 77 fermion 
[ 3 ]. In particular, time-reversal could be either Kramers 
(T^ = —1) or non-Kramers (T^ = !)• 

In 3D, band structures of Kramers fermions with T 
symmetry are classified by an integer index m V which 


counts the number of Majorana cones on the surface. It 
was realized [38114(1] that in the presence of interactions 
the state with tz = 16 is trivial, while that with tz = 8 
is equivalent to a topological paramagnet. More specif¬ 
ically, for 7 z = 8 the surface state with four Dirac cones 
(eight Majorana cones) can be gapped without breaking 
any symmetry via strong interactions, and the result¬ 
ing gapped surface state must have intrinsic topologi¬ 
cal order. The simplest such topological order is a Z 2 
gauge theory in which the e and m particles are bosons, 
but transform under time-reversal as Kramers doublets 
(T^ = — I). Therefore we can put the slave fermions into 
a band with z/ = 8 , and let the gauge fields confine the 
fermions (either automatically through an SU(2) gauge 
field or at strong coupling through a U(I) gauge field). 
Crucially, the topological quasi-particles (e and m) on 
the surface do not carry the gauge charge, and they sur¬ 
vive on the surface as deconfined objects. The resulting 
phases are therefore confined paramagnets with nontriv¬ 
ial surface states protected by time-reversal symmetry. 

Non-Kramers fermions, by contrast, cannot host non¬ 
trivial band structure with time-reversal symmetry alone. 
However, if spin- 5'2 conservation is present, the band 
structures can again be assigned an integer topological 
invariant v' m which is the number of Dirac cones on 
the surface (or half the number of Majorana cones). It 
is known [551 - HD] that with interactions the state with 
i/' = 8 is trivial, while that with v' = A \s equivalent to 
a topological paramagnet. We can then put the slave- 
fermions into a band with v' = A and let the gauge fields 
confine the fermions, which produces a topological para¬ 
magnet with time-reversal and spin -52 conservation. 

In both cases we need to put the slave fermions into 
band structures with four Dirac cones on the surface. 
Band structures with two Dirac cones {v = 4) have 
been studied on the cubic im and diamond [ 12 ] lattices. 
Therefore we can obtain the desired structure simply by 
putting the partons into two copies of the v = A band. 
This can be easily done by taking advantage of the two 
orbitals in Eq. [^ making the topological paramagnets 
very natural in spin -1 systems. 

In the next section we outline a similar construction 
for the one-dimensional Haldane chain, by confining slave 
fermions which form four copies of the Kitaev chain. This 
illustrates the essential idea of our constructions in a sim¬ 
pler and more familiar context. 


A. Parton construction for Haldane/AKLT chain 

The Haldane phase is an SPT phase with gapless 
boundary degrees of freedom that are protected by time 
reversal. As a warm-up exercise, we outline how this 
phase can be constructed from a topological supercon¬ 
ductor of slave fermions. This illustrates some features 
we will meet again in 3D. A different parton construction 
for the Haldane phase was considered in Ref. [421 

The fermions are taken to be non-Kramers (T^ = !)■ 
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In ID, superconducting band structures for free non- 
Kramers fermions are labelled by a Z-valued index m, 
V, which is the number of protected Majorana zero modes 
at the boundary. The state with a given v can be viewed 
as V copies of Kitaev’s p-wave superconducting chain 
[44] , Interactions reduce this classification to Zg, i.e. the 
iz = 8 state becomes trivial |^. Further, the state with 
iz = 4 is topologically equivalent to the Haldane chain, 
modulo the presence of gapped fermions in a trivial band. 

Here we therefore put the slave fermions into four 
copies of the Kitaev bandstructure, in an SU(2)- 
symmetric manner. Gauge fluctuations (or Gutzwiller 
projection) will then remove the unwanted degrees of 
freedom, leaving a topological paramagnet in the Hal¬ 
dane phase. 

Starting with an antiferromagnetic spin-1 chain. 


% = JY,S^.S,+^ + ..., (19) 

i 

we represent the spins with slave fermions as in Eq. 
or equivalently Eq. The valence bond picture of the 
AKLT state suggests using a mean-field Hamiltonian for 
the partons with hopping t and spin-singlet, orbital- 
singlet pairing A, 


Hmf = “X! fih+i 


+ h.c. 


A 




- h.c. 


In terms of the Majoranas, this is 


Hmf =M = tpV + iAp^a^ry. (20) 

i 


We first consider this as a free fermion problem, then 
include the gauge fluctuations. 

Eor simplicity take A = t, which makes the terms in 
Hmf for different links commute. The Hamiltonian is 
simply four copies of the Kitaev chain, as can be seen 
immediately by going to a basis where cr^r^ is diagonal. 
To be more explicit, it is useful to define the matrix 


X = 


( 21 ) 


since M = tp^(l -|- p^ayry) = {2tpy)'P+, etc. So we may 
rewrite Hmf as 

Hmf = 

i 

Taking open boundary conditions, and denoting the left¬ 
most site of the chain by L, we see that the four modes 
in do not appear in the Hamiltonian. 

These four Majoranas correspond to two complex 
fermion modes that can be occupied or unoccupied, i.e 
to a degenerate four-dimensional boundary Hilbert space. 
At the level of free fermions, this degeneracy is protected 
by time reversal symmetry T, under which is invari¬ 
ant (since by dehnition Xp^^^ = 

Once we go beyond mean field theory, the fermions 
are coupled to confining gauge fluctuations. We will see 
below that two of the four boundary states are not gauge 
invariant — i.e. they can be thought of as having an 
unscreened gauge charge sitting at the end of the chain. 
Confinement removes these states from the low energy 
Hilbert space, leaving a single boundary spin-1/2 whose 
gaplessness is protected by time reversal. 

Hmf treats spin and orbital degrees of freedom sym¬ 
metrically, and preserves SU(2)spin x SU(2)orbitai sym¬ 
metry. The four boundary states can be labeled by the 
occupation numbers of two complex fermions cip. Since 
the partons transform as doublets under each SU(2), the 
fermions cip should also form doublets under each SU(2). 
In an appropriate basis the transformations are 

SU( 2 )spin : (01,02)"^ ^ Z^s(ci,C 2 )^, 

SU(2)orbitai : (ci,4f ^ Ho(ci,4f. (26) 

where are SU(2) matrices. It follows that states 
which are singlets under SU( 2 )spin are doublets under 
SU(2)orbitai and vice versa. We denote the spin doublet 
It), li) and the orbital doublet | 1 ), | 2 ). The spin operator 
for the boundary spin -1 can be split into contributions 
from the dangling boundary modes p^^ and from p^ 

4 = 4 +^ + si-\ with 


Firstly, we use this to dehne the action of time reversal 
T on the fermions: 


^(±) ^ 1 (±)Tg^(±)^ s = (pV^,a^pV^). (27) 

8 


r: p-^Xp. (22) 

This definition ensures that the spin changes sign under 
T and that Hmf is invariant. The fermions are non- 
Kramers (T^ = 1 on 77 ). 

Secondly let us define matrices that project onto a 
given value of X, and corresponding fermion modes: 

r± = ^{i±x), 4^^=v±p. (23) 

In an appropriate basis, ryi+i has four nonzero compo¬ 
nents. Next, note that 

(24) 


We can make a similar splitting for the orbital spin T, 
which is related to S by swapping the as for rs. We 
denote the matrices appearing in T by D: 

f(±) ^ 1 (±)T^^(±)^ ^ ^ (pyT^,Ty,pyT^). (28) 

8 


^ Any quadratic term (where A is real antisymmet¬ 

ric) is forbidden as it is odd under T. However in the presence 
of interactions a four fermion term 71727374 — where ji are the 
components of in some basis — is allowed by time reversal, 
and lifts the boundary degeneracy to a single doublet as in the 
Haldane chain | 45 | . 


M = V-MV. 
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The pairs (|t) , li)) and (|1), |2)) are both Kramers dou¬ 
blets, since the spin and orbital operators for the bound¬ 
ary modes, and , change sign under T. This can 
also be checked explicitly by considering the transforma¬ 
tion of the boundary states (labeled by fermion occupa¬ 
tion numbers) under T, with the fermions transforming 

as T : Cy2 —t c| 2- 

Now we consider the effect of gauge fluctuations or 
Gutzwiller projection. We have listed the generators for 
the Sp(4) gauge group in Eq. 18 However, some gauge 
generators are Higgsed in the above mean field state. In 
general, to determine the unbroken gauge group we must 
examine Wilson loops of the form W = 
where 7 ?mf = Thij El- The unbroken gauge gen¬ 

erators are those that commute with the Wilson loops. 
Here, the only nontrivial Wilson loop is the matrix X 
defined in Eq. This leaves a subset of six unbroken 
generators, which may be written in terms of the matri¬ 
ces H appearing in the orbital spin (Eq. 281: 


riD = {H,^H}. (29) 

Taking linear combinations, we can use insteacQ 

Tiu = {v+ftv+,r-(ir-}. (30) 


It is interesting to consider inversion symmetry here. 
In the free fermion problem, v —> —v under inversion, 
so that a nonzero value of v can only be realised with 
a Hamiltonian which breaks inversion symmetry. With 
interactions, ~ v + 8, suggesting that u = A can be 
realised in inversion-symmetric interacting system [47j . 
The present example is a nice realisation of this. The 
mean field Hamiltonian Hmf appears to break inversion 
symmetry. However, the symmetry can be restored by 
combining it with a gauge rotation. So the projected 
wavefunction is actually inversion symmetric. 

We now move on to 3D states. 


B. Cubic lattice 

Making use of the cubic band structure studied in 
Ref. im we construct an SU(2) gauge theory which con¬ 
fines to a topological paramagnet. We choose the mean 
field Hamiltonian 

Hmf =^Uj'nI ( 31 ) 

{ij) {{ij)) 

(((y»> 


We denote the unbroken gauge group SU(2)^^j]; 

SU(2)^”^ 


orbital 


orbital' 


To make the Hamiltonian in Eq. (20) a reasonable 
ansatz, we must check that the Sp(4) gauge charges are 
all zero on average: (T^) = 0 for all i. Fortunately the 
unbroken gauge symmetry Tid guarantees this. 

The boundary modes involve only , so are invariant 
under SU(2)^~|]j^j^j. However, |1) and |2) are not invariant 
under «S'{7(2)^(('j]jjg^j. Therefore after confinement only the 

doublet It), It) survives, with corresponding spin 
This is the boundary spin-1/2 of the Haldane phase. 

In this ID example we can confirm explicitly that 
Gutzwiller-projecting the mean-field wavefunction gives 
the desired SPT phase. In fact the Gutzwiller-projected 
state for A = t, denoted 14'spin), is precisely the AKLT 
state. To see this we adopt a trick from Ref. 031 Using 
the fact that the terms in Hmf commute, we can check 
that l^^spin) has zero amplitude for a pair of adjacent sites 
to be in a spin-two state, l^/gpin) is therefore the ground 
state of the AKLT Hamiltonian, since this can be writ¬ 
ten as a sum of projectors onto the spin-two subspace for 
each linkEl 


where the nearest-neighbor hopping tij gives a tt-Aux on 
every square plaquette, the body-diagonal pairing Xij fol¬ 
lows the pattern studied in Ref. 011 and the next-nearest- 
neighbor pairing Xij is a small perturbation introduced 
to reduce the gauge group to SU(2) and is not responsible 
for the gap or the band topology. 

To determine the unbroken gauge group, we examine 
the Wilson loops as above. The fundamental nontrivial 
ones are proportional to and The unbroken 

gauge group is generated by those of the Sp(4) generators 
that commute with the Wilson loops. It is then straight¬ 
forward to see that the unbroken gauge group is an SU(2) 
generated by 


rcubic = {pV^r",r^pW^}. (32) 

One can choose to implement time-reversal T as T : 
r] —)• ip^a^r], and it is straightforward to see that T : 
Hmf ^ Hmf j S ^ and r^ubic ^ Tcubic- The band 
structure in Eq. (31) preserves time-reversal symmetry, 
and the SU(2) gauge rotation commutes with T. Notice 
also that = — 1 on the ry fermions. 


^ To be more precise, the two types of generators in Eq. |29| c or- 
respond to elements of the invariant gauge group dGG) I3l at 
different momenta, k = 0 and k = n, so when we take linear 
combinations the two types of generators in Eq. |30| alternate on 
even and odd sites. This is not crucial here. Another subtlety is 
that the IGG is enlarged at the special point A = f. 

® This correspondence with the AKLT state is less obvious if 


we simply Gutzwiller—project the BCS ground state of Hmf- 
The ground state of the Kitaev chain involves a long-range 
Cooper pair wavefunction, C{r) = {L — 2r)IL so in the 

present case j^'spin) is obtained by acting on the vacuum with 
an exponented sum of long-range singlet creation operators, 
exp(I]i l]r >0 UR)////((/), and projecting. The AKLT 
state may of course be written using only short-range singlet 
creation operators, |AKLT} oc V Yliifi fj^i) |vac}. By the 
previous argument, the two states must be equivalent. 
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We must check that the Sp(4) gauge charges are all 
zero on average, (Fi) = 0. The unbroken gauge sym¬ 
metry Fcubic guarantees that (r^) = 0 for all i except 
for Fs = . Furthermore, time-reversal invariance T 

guarantees that (F5) = 0. Hence the condition is indeed 
satisfied for any i. 

To determine the band topology, it is sufficient to 
consider the Hamiltonian with only the nearest- 

neighbor and body-diagonal terms in Eq. 
fermions with different orbital indices are decoupled and 
form two identical bands. Each band is the same as that 
studied in Ref. SH with V = A (two Dirac cones on the 
surface). So the band has :/ = 8 in total (four Dirac 
cones). So Eq. (31) indeed gives rise to a topological 
paramagnet. 

In order to understand the role played by spin- 
rotation symmetry, we examine the surface state in 
more detail. We start from the surface Dirac theory 
with SU(2)gauge X SU(2)spin X T symmetry, with four 
Dirac cones in total: 


H = 'ilj\pxPx+PyPz)‘^To‘^aoi>, (33) 
with time-reversal 


T Ip ^ ipy®TQ® O’oV’^ : 

(34) 

gauge SU(2) 


Ug : Ip ^ p.o®Ug® (jQip, 

(35) 

and spin SU(2) 


Us : Ip ^ po ® tq ® Usip. 

(36) 


We have denoted the SU(2)gauge Pauli matrices by r, but 
they should not be confused with the Pauli matrices for 
the orbital spin. 

Next we will consider driving this surface theory into 
a Z2 topologically ordered state by first introducing an 
order parameter A which gaps out the Dirac fermions, 
but breaks time reversal symmetry, and then restoring 
time-reversal symmetry by proliferating double vortices 
in A. The single vortex remains gapped, and gives rise to 
anyonic surface excitations with nontrivial time reversal 
properties. 

To analyse the symmetry properties it is useful to con¬ 
sider the auxiliary U(l)a transformation 

Uai9) : ^ ^ e*V (37) 


(which is an emergent symmetry of Eq. 33 but not a 
microscopic symmetry). The gap term of interest is 


Ha = lA ^l}py '^Ty ^ ffyip F h.c. (38) 

This is invariant under the SU(2)gauge x SU(2)spin sym¬ 
metry. It is not invariant under time reversal T or under 
U(l)a separately, but it is invariant under the modihed 


time-reversal transformation 7” =Ua{TT/2)T. Notice that 
= 1 on the parton fermions ip, in contrast to the orig¬ 
inal T under which they are Kramers. 

As shown in Refs. m iMi and EOl the fundamental 
vortex in A transforms projectively under T, i.e. = 
— 1. We now examine the SU(2)gauge x SU(2)spin spins 
carried by the vortex. A key point is that there are four 
Majorana zero modes trapped in the vortex core. One 
can label the internal Hilbert space with two complex 
fermions ci^2- Since both SU(2) groups are preserved in 
the intermediate gapped phase and the partons transform 
as doublets under both SU(2), the two complex fermions 
ci^2 should also be doublets under both SU(2). In an 
appropriate basis the transformations are 

Ug :(Cl,C2)^ Ug{Cl,C2)^ , 

Us :{ci,c\Y -^Us{ci,c\)'^. (39) 

It follows that states which are singlets under SU(2)gauge 
are doublets under SU(2)spin and vice versa. Specifically, 
there are two distinct kinds of vortices, labeled by the 
fermion parity (—l)'=i'=i+4'=2. both have = —1, but 
one transforms as (0,1/2) under SU(2)gauge x SU(2)spin, 
and the other as (1/2, 0). 

We now restore time-reversal symmetry by condens¬ 
ing double-vortices that transform trivially under both 
SU(2)gauge X SU(2)spin and T, giving Z 2 topological or¬ 
der on the surface [48l |49]. Single vortices with even 
and odd fermion parity yield mutual semions which 
we denote e and m respectively. Both are Kramers 
bosons (T^ = —1), and e transforms as (0,1/2) under 
SU(2)gauge X SU(2)spin while m transforms as (1/2,0). 
Their bound state, e, is non-Kramers, fermionic, and 
transforms as (1/2,1/2). 

So far, our treatment of the surface has neglected the 
confining gauge heldj^ When we take it into account, 
only excitations that are neutral under SU(2)gauge sur¬ 
vive. In addition to e, these include bound states m = ipe 
and e = iprh got by attaching a ip fermion to m and e. 
This shifts the self-statistics, so m is bosonic while e is 
fermionic (all three particles are mutual semions). Since 
e is the bound state of e and m (and its properties follow 
from this) we do not discuss it further. Note that m = ipe 
is Kramers since ip is. 

The upshot is that the surface topological order sur¬ 
viving after ‘gauge neutralization’ has an e particle that 
is Kramers and spin-doublet, and an m particle that 
is Kramers but spin-singlet. Since both e and m are 


The fermions in the bulk are confined, giving a non-fractionalized 
bulk state. It is known m that the SU(2) gauge theory has a 
0-term &t 9 = {u'I2)7 t. For u' = 4, as here, we have 6 = 27r 
which has the same physics as at 0 = 0. The confined state then 
can preserve time reversal symmetry. In contrast if = 2, we 
will have 0 = ir and the resulting confined phase of the SU{2) 
gauge theory must be non-trivial in some way. It either breaks 
time reversal or becomes a quantum spin liquid with long range 
entanglement. 
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Kramers bosons, this state is indeed the eTmT phase, 
like the wavefunction discussed in Sec. HI 

However if spin-rotation symmetry is preserved, a finer 
classification is possible, under which the present state is 
dubbed eCTmT, where the ‘C’ indicates that e is a spin 
doublet [To] Q This finer classification emphasises a differ¬ 
ence between the eTmT state constructed here, in which 
e is a spin doublet and m is not, and that constructed in 
Sec. U where both e and m are spin doublets. 

Like the ID example of the previous section, the cu¬ 
bic lattice construction violates inversion symmetry at 
the free fermion level (this is inevitable if v is nontrivial 
M) but the resulting spin state is inversion symmet¬ 
ric as a result of gauge invariance. Here, the hopping 
term in Hmf is invariant under inversion, while for an 
appropriate choice of x' th® pairing terms change sign 
under inversion. Therefore inversion can be restored by 
combining it with the gauge transformation / —> i/, i.e. 
r] —)■ ip'^rj. (With the arrow conventions of Sec. IB the 
fluctuating AKLT state is also inversion symmetric.) 


C. Diamond lattice 


Next we consider parton theories on the diamond lat¬ 
tice, making use of the band structure of Ref.jHj First we 
construct a theory with an SU(2) gauge held which natu¬ 
rally conhnes (Sec. jlH Cl ). The resulting state is a topo¬ 
logical paramagnet which requires both time-reversal and 
AK-spin rotation symmetry to be protected. Then in 


Sec. HI C 2 we construct a U(l) gauge theory, which con¬ 


hnes at strong coupling. The conhned state is a topolog¬ 
ical paramagnet which only requires time-reversal sym¬ 
metry. 


1. Topological XY paramagnet from SU(2) gauge theory 
The mean held Hamiltonian is 

Hmf = pypj+ ^ t'ijvTp^Vj+ 

(ij) ((b» ((b» 

(40) 

where the nearest-neighbor hopping t is isotropic, while 
the next-nearest-neighbor hopping fb and pairing 
follow the patterns discussed in Ref. Notice that the 
pairing term is a singlet in orbital space, but is a triplet in 
spin space. Hence the spin rotation symmetry is reduced 
from SO(3) down to 0(2) rotations about the Sy axis, 
corresponding to XY anisotropy in the spin model. 

We again calculate the nontrivial Wison loops: the 
simplest nontrivial ones consist of three links and are 


^ The eCTmT state is topologically equivalent to the combination 
of a generic eTmT state and the state eCmT\ the analogue of 
the latter for U(l) spin symmetry is discussed in Sec. IIII C 1| 


proportional to p^ and p^r^. The unbroken gauge group 
is generated by 

^diamond = {p*'T^T^p«T"}. (41) 


These are precisely the orbital SU(2) generators H. 

One can implement time-reversal symmetry T as ry —>■ 
p^tT^r^Ty, under which t] is non-Kramers (T^ = 1) and 
S ^ —S and of course Hmf Hmf- 

As above we must check that the Sp(4) gauge charges 
are all zero on average: (Ti) = 0. The unbroken gauge 
symmetry Tdiamond guarantees that (T^) = 0 for all i 
except for Ti = p^, which is nothing but the total fermion 
occupation number (minus two). Fortunately the mean 
field Hamiltonian Eq. (401 has a special lattice symmetry 
m that sets (Fi) =0. 

To determine the topology of the mean field band 
structure, it is convenient to consider the modified time- 
reversal symmetry T' : p -T ip’^T^rj (with = —1), 
which is the combination of time-reversal and spin rota¬ 
tion ioy. Fermions with different physical spins (p^ and 
rji) do not mix under the modified time-reversal. Fur¬ 
thermore, they are decoupled in the mean field Hamil- 
tionian Hmf and form two copies of an identical band. 
Therefore the topological index v' is defined for each 
band separately. Now each band is identical to that stud¬ 
ied in Ref. HU with i/' = 4. The total band therefore has 
ly' = 8, with four Dirac cones in total on the surface. 

We now consider the surface Dirac theory 
with SU(2)gauge X U(l)spin X T Symmetry, with 
four Dirac cones in total: 


H = -iPHpxTx + PyTz) ® T-q (g) CToV', 

(42) 

with modihed time-reversal 


T' ■■ tp ^ ipy ® Tq ® croV’^ 

(43) 

gauge SU(2) 


Ug Ip ^ p,o®lylg® O’oV': 

(44) 

and spin U(l) 


C/«(0) : V'^ e*V. 

(45) 

The actual time-reversal is T = U/‘^)T''■ 


T ■.'p^pLy®TQ® (TQlp^ 

(46) 

Now consider the gap term 


H^ = lA ^PPy ®Ty® a yip -|- k.C., 

(47) 

which preserves both SU(2)gauge and T, but breaks 
U(l)spin- To restore the U(l)spin symmetry and preserve 
the gap, we need to proliferate vortices in the order pa¬ 
rameter field A. It was shown in Ref. IT^ [39l and HO] that 
the fundamental vortices have = —1, so condensing 


double-vortices gives a Z 2 gauge theory, with e being 
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Kramers, m being Kramers and SU(2)gauge-doublet, and 
e being non-Kramers and SU(2)gauge-doublet. We can 
then gauge-neutralize the particles by binding ip fermions 
to rh and e. The neutralized theory then has e being 
Kramers, and m = etp being non-Kramers (recall that 
= 1 on ip) but carrying spin-1/2 under U(l)spin due 
to the S„ spin carried by ib. This state is dubbed eCmT 

inRef.Iini 

The fermions will be confined once the fluctuation of 
the SU(2) gauge field is introduced, and we obtain a 
non-fractionalized bulk state. On the surface, the eCmT 
topological order survives the confinement, since all the 
non-trivial quasi-particles in the theory are gauge-neutral 
and are hence decoupled from the gauge field. We have 
thus obtained the eCmT topological paramagnet. 

As a side note, if the spin-1 operators are pseudo-spins 
such that T : Sy, S^} —>■ {^a,, —S'y, 52}, then the 

modified time-reversal T' '■ rj —>■ ip^r’^rj (with = ~1) 
could represent the physical time-reversal symmetry. In 
this case we obtain a topological paramagnet that re¬ 
quires time-reversal only, as will be shown in Sec |III C 2| 


2. Stable U{1) quantum spin liquids and topological 
paramagnets 

The parton construction of course also gives access to 
stable quantum spin liquid phases. Of particular interest 
to us is a time reversal symmetric U(l) quantum spin liq¬ 
uid phase on the diamond lattice. For greater generality 
we allow for full SU(2) spin symmetry. As usual such a 
phase has a gapless emergent photon. In addition it has a 
gapped fermionic spin-1/2 Kramers doublet spinon which 
has internal ‘electric’ charg^and a gapped bosonic spin- 
0 magnetic monopole that transforms to an antimonopole 
under time reversal. We will give the spinons the band 
structure of a topological superconductor (as in previ¬ 
ous sections). The resulting quantum spin liquid phase 
then inherits the non-trivial surface states of the topolog¬ 
ical superconductor. The relevance to the present paper 
comes from asking about the confined phase that results 
when the magnetic monopole is condensed. We show be¬ 
low that this is the eCTmT topological paramagnet. 

SPT phases in 3D have been accessed previously 
through confinement of emergent U(l) gauge fields [35] . 
However in these previous studies the confinement was 
achieved in a highly non-trivial way involving the con¬ 
densation of dyons (bound states of magnetic and elec¬ 
tric charges). The novel aspect of our construction is 
that the confinement is achieved directly by simply con¬ 
densing the magnetic monopole, which will result from 
the usual dynamics of the gauge field at strong coupling. 


This ‘electric’ charge couples to the emergent photon in this spin 
liquid, and not to physical external electromagnetic fields. 


Consider the following mean field ansatz: 

Hmf =^tr]Jpyr]j + ^ -I- A.jTifp^aypj 

(b> <<b)) (<b>) 

+ (48) 

i {ij) A 


where the nearest-neighbor hopping t and on-site pair¬ 
ing A' are uniform and isotropic, while the next-nearest- 
neighbor hopping fb and pairing Ay follow the patterns 
discussed in Ref. [421 Note that the first two terms are 
the same as in Eq. ( [40| ), and the third is got by exchang¬ 
ing the role of orbital and physical spin. Contrary to 
Eq. (40), the pairing term A is a singlet in physical spin 
and a triplet in orbital space, so the full spin rotation 
symmetry is preserved. The nearest-neighbor antisym¬ 
metric hopping term t" is introduced to reduce the gauge 
symmetry, and does not affect the other arguments in this 
section as long as it is kept small. 

The simplest nontrivial Wilson loops are proportional 
to p^, p^a^ and p^T^. The resulting unbroken gauge 
group is a U(l) generated by r^. 

We implement time-reversal symmetry T through rj —>■ 
ip^a'^r] (which has = —1). It is straightforward to 
check that S —>■ —S and Hmf —> Hmf under the cho¬ 
sen time-reversal symmetry. Moreover, the U(l) gauge 
charge is also odd under T, which allows for topolog¬ 
ically non-trivial band structures for the partons. 

We now check that (T^) = 0. The unbroken U(l) gauge 
symmetry and time-reversal guarantee that (T^) =0 for 
all i except for p^ and p^a'^, which are nothing but the 
total fermion occupation number (minus two) and the 
real part of the on-site pairing. The lattice symmetry 
m again sets (p^) = 0. For the on-site pairing ampli¬ 
tude, there is no symmetry to set it to zero automatically. 
We must therefore adjust the on-site pairing term A' in 
Eq. ( [4^ to make it zero on average [52] . 

To determine the topology of the mean field band 
structure, notice that fermions with different orbital in¬ 
dices (r indices) do not mix under time-reversal T : p —>■ 
ip^a^rj. They are also decoupled in the mean field Hamil- 
tionian Hmf, forming two copies of an identical band. 
Therefore the topological index i/' is defined for each 
band separately. Now each band is almost identical to 
that studied in Ref. |42l with 1 /' = 4. The total band 
therefore has = 8, with four Dirac cones in total on 
the surface. 

We now consider fluctuations of the U(l) gauge field. 
In the weak coupling regime the gauge theory is decon- 
fined, and we have a stable U(l) quantum spin liquid 
phase. The spinon band structure has time reversal pro¬ 
tected surface states that provide a distinction between 
this spin liquid and more conventional ones. For a com¬ 
pact U(l) gauge theory, there are always gapped mag¬ 
netic monopole excitations in the theory. In Ref. [T3| and 
1391 it was shown that for the spinon band structure we 
have here, this magnetic monopole is a spin-0 boson that 






14 


simply transforms into an antimonopole under time re¬ 
versal. 

As the gauge coupling strength increases, the 
monopole mass gap decreases and eventually becomes 
zero. The monopoles will then condense and confine the 
gauge theory. The trivial symmetry properties of the 
monopole implies that this condensate does not break 
T or the physical spin SU(2) (if present). The con¬ 
fined state is thus a non-fractionalized symmetry pre¬ 
serving paramagnet. To determine which SPT phase the 
paramagnet belongs to, we need to examine the surface 
state in more detail. The argument is largely parallel to 
that in Sec. IIIB with the simple modification that the 
SU(2) gauge symmetry discussed in Sec. IIIB is reduced 
to U(l). The conclusion remains the same: the param¬ 
agnet is the nontrivial SPT dubbed eCTmT in Ref. ITOl 
The representative surface state is a gapped Z 2 topologi¬ 
cal order, with e being Kramers and spin-doublet, and m 
Kramers but spin-singlet. (If the spin-rotation symmetry 
is broken, this becomes a generic eTmT state.) 


D. Spin wavefunctions 


The parton constructions suggest spin wave functions 
that may be useful as variational states in future work 
on specific microscopic models. Following the standard 
procedure [3] we construct a spin wave function from the 
mean-field fermion wave function I^^mf) by projecting 
onto the subspace obeying the constraints J2aa flafaa = 
2 and Y.aba fL^abha = 0: 

l^-spin) =iP|«'MF). (49) 


Such a projection is expected to roughly mimic the ef¬ 
fect of gau ge f luctuati ons. For the states constructed in 
Sec. IIIB and IIIC 1[ the SU(2) gauge fluctuations au¬ 
tomatically confine the states. We therefore expect the 
projected wave functions to repr esent th e confined spin 
SPT states. For the state in Sec. IIIC2 the U(l) gauge 
field is deconfined at weak coupling, and confines to an 
SPT state at strong coupling. So it is not clear a priori 
whether the projected wave function will give the U{1) 
quantum spin liquid state or the confined SPT state. 

These spin wave functions are alternate possibilities to 
the loop gas wave functions described in the first part 
of the paper. While the loop gas wave functions are 
physically appealing they are likely not very tractable 
numerically due to the linking signs. The parton wave 
functions, on the other hand, may be studied through 
variational Monte Carlo calculations though the physical 
connection to SPT physics is less directly obvious. This 
situation is similar to existing descriptions of quantum 
spin liquid phases through either loop gases (string-nets) 
or through partons which each have their advantages and 
disadvantages. 

For the topological paramagnets, at present we do not 
have a direct connection between the parton and loop gas 
wavefunctions. Establishing such a connection is a target 


for future work, and will confirm the general correctness 
of the projected wave functions as faithfully capturing 
the state accessed through the parton description. 


IV. DISCUSSION: TOWARDS MODELS AND 
MATERIALS 

We have emphasized that frustrated spin-1 magnets in 
3D may be fruitful in the search for spin SPT phases. 

In the ongoing search for quantum paramagnetism in 
frustrated systems, the bulk of the attention has focused 
on spin-1/2 systems. This is guided by the intuition that 
increasing the spin only leads to more ‘classical’ physics 
and hence to a greater tendency to order. Caution how¬ 
ever is required in taking this intuition too seriously. In 
one dimension the spin-1/2 chain is almost antiferromag- 
netically ordered (power law correlations) while the spin- 
1 chain is a good paramagnet with a spin gap. This 
has the following amusing consequence. Consider a two- 
dimensional rectangular lattice with nearest neighbor an¬ 
tiferromagnetic interactions: 

Hrect — J\\ ^ ^ Sr ■ Sr+x. + J± E Sr • Sr^y (50) 

r r 

For J|| = J_L the model is antiferromagnetically ordered 
for all spin S. When ^ is decreased from 1 the spin-1/2 
model stays ordered unless J_l = 0. The spin-1 model 
on the other hand becomes a spin gapped paramagnet 
below a non-zero critical value of • So there is a range 
of parameters in this 2D model where the spin-1 system 
is a quantum paramagnet although the spin-1/2 system 
has long range Neel order. 

There are some interesting examples of frustrated spin- 
1 magnets — most notably NiGa 2 S 4 and Ba 3 NiSb 209 , in 
both of which the spin-1 Ni ion forms a triangular lattice 
[531154] . Apart from new and interesting kinds of quan¬ 
tum spin liquids, spin-1 magnets may also harbor novel 
broken symmetry states (such as spin nematics [55] 1 more 
naturally than their spin-1/2 counterparts. To this we 
add the SPT phase discussed in this paper as a possible 
fate for a frustrated 3D spin-1 magnet. 

Our results suggest a route to guessing possible micro¬ 
scopic models that might harbor an SPT phase. Start¬ 
ing from the parton mean field Hamiltonian we can write 
down a lattice gauge theory that captures fluctuations. 
A strong coupling expansion of this lattice gauge the¬ 
ory will result in a spin Hamiltonian which may then be 
in the same phase as the same lattice gauge theory at 
weaker coupling. Such an approach has previously been 
successfully used to write down lattice models for various 
spin liquid phases. Given that we are interested here in 
confined phases we may be cautiously optimistic that a 
similar approach has an even better chance of resulting 
in spin models for the SPT phases. As an application 
let us consider the diamond lattice parton construction. 
With full SU(2) spin symmetry, the mean field state of 
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Section IIIC 2 suggests (at leading order of the strong 


coupling expansion in the resulting U(l) gauge theory) 
an interesting frustrated spin-1 model: the “J1-J2' 
tiferromagnet on the diamond lattic^ 


an- 


H — Ji Sr ■ Sr' + J2 ^ Sr ■ Sr 


( 51 ) 


{rr') 


{{rr')) 


The next-nearest neighbour coupling J2 introduces frus¬ 
tration. Indeed classically once J2 > ^ there are an 
infinite number of degenerate ground states | 56 | that are 
not related by global spin rotation. For large spin, it 
has been argued that the ground state is magnetically 
ordered as a result of quantum order by disorder m- 
The ground state for S' = 1 (or S' = 1 / 2 ) is not known. 
The SPT paramagnet discussed in this paper is a candi¬ 
date. The various descriptions we have provided should 
be a useful guide in future numerical studies should a 
paramagnetic ground state be found for this model. 

It is interesting to note that since the diamond lat¬ 
tice is 4 -fold coordinated classical 2 -sublattice Neel or¬ 
der is likely to be more easily destabilized by frus¬ 
tration/quantum fluctuations than in the cubic lattice. 
Thus the J1-J2 diamond magnet for low spin [S = 1/2 
or 1) may be an excellent candidate to find an interesting 
quantum paramagnetic ground state. 

The frustrated diamond lattice model appears to de¬ 
scribe well [ 33 ] the physics of the spinel oxide materials 
MnAl2 04 and C0AI2O4 [ 58 ] which belong to a general 
family of materials of the form AB2O4. The A site forms 


the diamond lattice and is magnetic. The Mn and Co 
compounds have S = ^ and S' = | respectively. In 
searching for a material that realizes the S = 1 model 
it is natural then to consider NiAl2 04. Here Ni is ex¬ 
pected to be in a d® Ni^“'' configuration and have spin-1. 
However the A site is tetrahedrally coordinated, and in 
the resulting crystal field, the Ni^"'' ion will have orbital 
degeneracy in addition to spin- 1 . Further spin-orbit cou¬ 
pling will split the resulting spin-orbital Hilbert space 
and the physics of the lattice will be determined by its 
competition with inter-site spin/orbital exchange. Thus 
NiAl204 will not be simply described by a spin-1 dia¬ 
mond lattice model. Work toward obtaining an appro¬ 
priate ‘spin-orbital’ model for NiAl204 is currently in 
progress | 59 | . It remains to be seen whether the pres¬ 
ence of the orbital degrees of freedom aids or hinders the 
formation of paramagnetic states. 

In any case we hope that these considerations motivate 
an experimental search for and study of frustrated spin-1 
magnets. 
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